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ABSTRACT 

We study the effect of vacuum polarization on the atmosphere structure and radi- 
ation spectra of neutron stars with surface magnetic fields B = 10^"* — 10^^ G, as 
appropriate for magnetars. Vacuum polarization modifies the dielectric property of 
the medium and gives rise to a resonance feature in the opacity; this feature is narrow 
and occurs at a photon energy that depends on the plasma density. Vacuum polariza- 
tion can also induce resonant conversion of photon modes via a mechanism analogous 
to the MSW mechanism for neutrino oscillation. We construct atmosphere models in 
radiative equilibrium with an effective temperature of a few x 10^ K by solving the 
full radiative transfer equations for both polarization modes in a fully ionized hydro- 
gen plasma. We discuss the subtleties in treating the vacuum polarization effects in 
the atmosphere models and present approximate solutions to the radiative transfer 
problem which bracket the true answer. We show from both analytic considerations 
and numerical calculations that vacuum polarization produces a broad depression in 
the X-ray flux at high energies (a few keV <E<a few tens of keV) as compared to 
models without vacuum polarization; this arises from the density dependence of the 
vacuum resonance feature and the large density gradient present in the atmosphere. 
Thus the vacuum polarization effect softens the high energy tail of the thermal spec- 
trum, although the atmospheric emission is still harder than the blackbody spectrum 
because of the non-grey opacities. We also show that the depression of continuum flux 
strongly suppresses the equivalent width of the ion cyclotron line and therefore makes 
the line more difficult to observe. 

Key words: magnetic fields - radiative transfer - stars: atmospheres - stars: magnetic 
fields - stars: neutron - X-rays: stars 



1 INTRODUCTION 

Thermal radiation from the surface of isolated neutron stars 
(NSs) can provide invaluable information on the physical 
properties and evolution of NSs. In the last few years, such 
radiation has been detected in four types of isolated NSs: 
(1) radio pulsars (see Becker 2000), (2) old, radio-quiet 
NSs (not associated with supernova remnants), some of 
which may be accreting from the interstellar medium (see 
Treves et al. 2000), (3) soft gamma-ray repeaters (SGRs) 
and anomalous X-ray pulsars (AXPs) (see Hurley 2000; Is- 
rael, Mereghetti, & Stella 2001; Mereghetti 2001; Perna et 
al. 2001), which form a potentially new class of NSs ("mag- 
netars" ) endowed with superstrong {B > 10^* G) magnetic 
fields (see Thompson & Duncan 1996; Thompson 2001), 
and (4) young, radio-quiet NSs in supernova remnants, 
which may include AXPs or SGRs (see Pavlov et al. 2001). 
The NS surface emission is mediated by the thin atmo- 



spheric layer (with scale height ~ 0.1 — 10 cm and density 
~ 0.1 - 100 g/cm^) that covers the stellar surface. There- 
fore, to properly interpret the observations of NS surface 
emission and to provide accurate constraints on the physi- 
cal properties of NSs, it is important to understand in detail 
the radiative properties of NS atmospheres in the presence 
of strong magnetic fields (see Pavlov et al. 1995 and Ho & 
Lai 2001 for more detailed references on observations and 
on previous works of NS atmosphere modeling). 

This paper is the second in a series where we system- 
atically investigate the atmosphere and spectra of strongly 
magnetized NSs. In Ho & Lai (2001, hereafter Paper I), we 
constructed self-consistent NS atmosphere models in radia- 



tive equilibrium with magnetic field B 



and Toff ~ 10" - 10^ K and assuming the atmosphere is 
composed of fully ionized hydrogen or helium. The radia- 
tive opacities include free-free absorption and scattering by 
both electrons and ions computed for the two photon po- 
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larization modes in the magnetized electron-ion plasma. It 
was found that, in general, the emergent thermal radiation 
exhibits significant deviation from blackbody, with harder 
spectra at high energies; the spectra can also show a broad 
feature {AE/Esi ~ 1) around the ion cyclotron resonance 
Esi = 0.63 (Z/yl)(B/10" G) keV, where Z and A are the 
atomic charge and atomic mass of the ion, respectively (see 
also Zane et al. 2001), and this feature is particularly pro- 
nounced when Ebz ^ SkTcft (however, as we show in this 
paper, vacuum polarization can significantly diminish the 
equivalent width of the ion cyclotron feature; see Sections ^ 
and||). 

In this paper, we study the effect of vacuum polariza- 
tion on the atmosphere structure and radiation spectra of 
strongly magnetized NSs. It is well-known that polarization 
of the vacuum due to virtual e+e" pairs becomes signifi- 
cant when B>Bq, where Bq = rule's/eh = 4.414 x 10^^ G 
is the magnetic field at which the electron cyclotron en- 
ergy hujBe = TieB/nicC equals mc(? ■ Vacuum polarization 
modifies the dielectric property of the medium and the po- 
larization of photon modes (e.g., Adler 1971; Tsai & Er- 
ber 1975; Gnedin, Pavlov, & Shibanov 1978; Heyl & Hern- 
quist 1997b), thereby altering the radiative scattering and 
absorption opacities (e.g., Meszaros & Ventura 1979; Pavlov 
& Gnedin 1984; see Meszaros 1992 for review). Of particular 
interest is the "vacuum resonance" phenomenon, which oc- 
curs when the effects of the vacuum and plasma on the linear 
polarization of the modes cancel each other, giving rise to 
"resonant" features in the radiative opacities (e.g., Pavlov 
& Shibanov 1979; Ventura, Nagel, & Meszaros 1979; Buhk 
& Miller 1997). At a given density p, the vacuum-induced 
resonance is located at the photon energy 

Sv ^ 1.02 ( ^ ^ (^i^) f{B) keV, (1.1) 



1 g cm" 



10^* G 



where Yc = Z/A is the electron fraction and ,f{B) is a slowly- 
varying function of B [f{B) = 1 for B ^ Bq and is of order 
a few for B ~ lO" - 10^'' G ; see Section H . It has been sug- 
gested that this vacuum resonance may manifest itself as 
absorption- like features in the spectra of X-ray pulsars and 
magnetars (e.g., Ventura, Nagel, & Meszaros 1979; Bulik & 
Miller 1997). However, as we show both analytically and nu- 
merically in this paper (see also Bezchastnov et al. 1996), 
since the atmosphere spans a wide range of densities and 
Ey depends on p, the spectral feature associated with Ey is 
spread out significantly and manifests as a broad "depres- 
sion" in the radiation spectrum. 

Furthermore, in a recent paper (Lai & Ho 2002a, here- 
after Paper II), we study the effect of the mode conver- 
sion that is associated with the vacuum resonance (see also 
Gnedin et al. 1978; Pavlov & Gnedin 1984): a photon prop- 
agating outward in the NS atmosphere can convert from one 
polarization mode into another as it traverses the resonant 
density py, which is given by 



pv « 0.96 



(1.2) 



(Across the resonance, the orientation of the polariza- 
tion ellipse rotates by 90°, although the helicity does not 
change.) This resonant mode conversion is analogous to the 
Mikheyev-Smirnov-Wolfenstein (MSW) effect for neutrino 
oscillations (e.g., Bahcall 1989; Haxton 1995) and is effec- 



tive for E> a. few keV (for which the propagation is adia- 
batic) and ineffective for E < a, few keV (the adiabatic con- 
dition breaks down). Because the two photon modes have 
vastly different opacities, this vacuum-induced mode con- 
version can significantly affect radiative transfer in magne- 
tar atmospheres (see Section 2.4). To properly account for 
this mode conversion effect in the atmosphere models, one 
must go beyond the modal description of the radiation field 
by formulating and solving the transfer equation in terms 
of the photon density matrix and including the effect of the 
refractive index. This is beyond the scope of this paper. In- 
stead, we consider two limiting cases in this paper: no mode 
conversion and complete mode conversion. These two limits 
correspond to different ways of labeling the photon modes, 
and the resulting model atmosphere spectra should approx- 
imately bracket the true spectra. 

There have been few previous works studying the effect 
of vacuum polarization on the radiation spectra of strongly 
magnetized NSs. Bezchastnov et al. (1996) reported NS at- 
mosphere models with B ~ 10^" - lO" G and Tofi « 10** K. 
This temperature is not appropriate for SGRs and AXPs, 
which have Tch ~ a few x lO" K (see, e.g., Mereghetti 2001; 
Perna et al. 2001). Although the vacuum-induced mode con- 
version was not discussed, the labeling of the photon modes 
in Bezchastnov et al. (1996) amounts to assuming complete 
mode conversion (see Section 2.4). Bulik & Miller (1997) 
considered how blackbody radiation with T ~ 10* K is mod- 
ified as it passes through a tenuous plasma where the opac- 
ity is affected by vacuum polarization, but no self-consistent 
atmosphere modeling, which is needed to determine the den- 
sity and temperature profiles, was attempted, and mode con- 
version was neglected. Ozel (2001) studied magnetar atmo- 
sphere models (with B > lO" G and Tes ~ a few x 10** K) 
which include the effect of vacuum polarization on the opaci- 
ties (see also Shibanov et al. 1992 for B ~ lO^'^ G atmosphere 
models which include a similar vacuum polarization effect on 
the opacities) . Besides neglecting resonant mode conversion, 
Ozel's work also neglects the ion effect and adopts a "satu- 
ration" scheme to smooth out the vacuum resonance feature 
in the opacity. Our analytical and numerical consideration 
of the vacuum resonance feature shows that such a "satu- 
ration" scheme cannot be justified (see Sections ^ and 6.3). 
Indeed, as we discuss in this paper, even when neglecting res- 
onant mode conversion, great care must be taken to properly 
handle the narrow and density-dependent vacuum-induced 
opacity feature in the radiative transfer (see Section |^ . 

In Section ^, we study the effect of vacuum polarization 
on the photon propagation modes in a magnetized electron- 
ion plasma and clarify the nature of the vacuum resonance 
and the associated mode conversion phenomenon. In Sec- 
tion ^, we discuss the change in the radiative opacities as 
a result of vacuum polarization. Section ^ shows the qual- 
itative effect of the vacuum resonance feature on the NS 
surface emission. The difficulties encountered in solving the 
radiative transfer equation caused by vacuum polarization is 
discussed in Section ^ We present atmosphere models and 
spectra for different magnetic field strengths in Section ^ 
Section ^ summarizes and discusses the implications of our 
results. Several "toy" atmosphere models that include opac- 
ities which mimic the vacuum polarization effects are dis- 
cussed in Appendix ^ 
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2 PHOTON POLARIZATION MODES IN A 
MAGNETIZED ELECTRON-ION PLASMA 
INCLUDING VACUUM POLARIZATION 

In this section, we outline the derivation of the photon polar- 
ization modes in a magnetized electron-ion plasma including 
the eflFect of vacuum polarization. Following the standard 
convention, the two normal modes are termed the extraor- 
dinary mode (X-mode, j = 1), which is mostly polarized 
perpendicular to the k-B plane, and the ordinary mode (O- 
mode, j = 2), which is mostly polarized parallel to the k-B 
plane, where k is the wave vector along the wave propagation 
direction and B is the external magnetic field (e.g., Meszaros 
1992). We shall see that this classification of the modes be- 
comes ambiguous near the vacuum res onan ce, where mode 
conversion can occur (Sections 2.3 and 2.4). 



2.1 Dielectric and Permeability Tensors 

For a cold electron-ion plasma, the plasma contribution to 
the dielectric tensor e'^' can be written, in the coordinate 
system x'y'z' with B along z' , as (e.g., Shafranov 1967) 



where 







ig 
e 

T] 



1 - Mui 



s 

Xs 



E 



/ 1/2 

sign(gs)Ma V. 

Xs Us 



•(1-Mi)(l-Me) 



1/2 

(1-Mi)(l-Me)' 



(2.1) 
(2.2) 

(2.3) 
(2.4) 



In equations ([2.2|)-(2.4), the sums run over each charged 
particle species s (electron and ion) in the plasma, and 
Us = uj%s/'^^ and Vs = i^jp^/o;^, where ujbs = \qs\B/{msc) 
is the cyclotron frequency and uips = {A-KUsq^ / msY^'^ is the 
plasma frequency of charged particle s. Damping of the par- 
ticle motion is accounted for by As = 1 + ivsjuj, where 
Vs is the d amp ing rate. In the second equalities of equa- 
tions (|2!^-(2.4), we have taken the M = Am^liZm^) = 
Ucjvi — [uc/uiY^'^ ^ 1 limit and assumed small damping 
[vs <^ u or \s —* 1). The relevant dimensionless parameters 
are 



2 



(2.5) 



where the electron cyclotron frequency uJse, the ion cy- 
clotron frequency ujEi, and the electron plasma frequency 
ui-pe are given by 



Eue = hiUBe = h— = 1.158 Bi4 MeV 
rricC 



ZeB 



EBi = hLOBi^ fi^P^ = 0.6305 Bi4 (^-jj keV 



AnipC 



(2.6) 
(2.7) 



Epe = hijpe = h 



t 28.71 ^^J 



/ 7\ 1/2 
f Z\ ' 1/2 



eV, 



\ rric 

respectively, B14 = B/lW^"^ G), and pi — p/(l g cm~^). 



Vacuum polarization contributes a correction to the di- 
electric tensor: 



Ae 



(v) 



(a - l)I-|-gBB, 



(2.9) 



where I is the unit tensor and B — B/B is the unit vector 
along B. The magnetic permeability tensor ^ also deviates 
from unity because of vacuum polarization; the magnetic 
field H„ of an electromagnetic wave is related to its mag- 
netic induction Sw by 
1 



H 



(al + mBB) ■ B^ 



(2.10) 



For hu) -C meC^, general expressions for the vacuum po- 
larization coefficients a, q, and m are given in Heyl & 
Hernquist (1997a). In the weak-field limit, B <C Bq — 
4.414 X 10^^ G, they are given by 



a = 
q = 

m — 
where 
Sv 



1 ~ 2Sv 

7<5v 

-45v, 



a^^2 
457r ' 



(2.11) 
(2.12) 
(2.13) 

(2.14) 



Qp = /he — 1/137 is the fine structure constant and 
b = B/Bq (Adler 1971). For B>Bq, we use the expan- 
sions given in Heyl & Hernquist (1997a,b; see also Tsai & 
Erber 1975) to find 



1 + ^ 
27r 



1.195 - I In 6 - i (0.8553 + In b) 



1 

'262 



'2^ 



(0.3070 + In 6) 



--6 + 1.272 - 
-0.70037^ 
1 + 1(0.1447 -ln6)-l 



(2.15) 

(2.16) 
(2.17) 



When I Ae<j' ] < 1 or 6 < Syr/ap (B < 5 x lO'^ G), the 
plasma and vacuum contributions to the dielectric tensor 
can be added linearly, i.e., e = e'''-' + Ae^^'. In the frame 
with B along z' , 
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(2.18) 










i . 










where a = a — 


1, e' = 


e + a, and rj 




+ a + g. The total 



magnetic permeability is given by equation (2.10) 



2.2 Equations for the Polarization Modes 

Usin g the electric displacement D — e ■ E and equa- 
tion ( ^.lOj ) in the Maxwell equations, we obtain the equation 
for plane waves with E oc e»(fc '— '^'): 

1 



+ n 



■J ^ij 



-{k X B)^{k X B) 



E, = 0, (2.19) 



where n = ck/uj is the refractive index and k = k/k. In the 
coordinate system xyz with k along the z-axis and B in the 
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x-z plane, such that k x B — — sin^sy (O b is the angle 
between k and B), the dielectric tensor [eq. (2.18)] is given 
by 

e' cos'^ 6's + 77' sin'^ ig cos 9b 

—ig cos 9b s' 
{e' — r]') sin 9b cos 9b ig sin 9b 

ie' ~ rj') sin9B cos 9b 



—ig sin 9b 
e' sin? 9b + ri' cos^ 9b 



The z-component of equation ( 2.19 ) gives 

Ez = —^zz i^zxEx + ezyEy) . 



Reinserting this back into equation (2.19) yields 
^ ' Ex 



Vyx Vvv ~ ''"^ J V 



= 0, 



where r = l + f = l + (m/a) sin^ 9b and 

^xx 



Vyy 



1 . _ N 1 ' ' 

y^zz^xx ^xz^zxj — ^ ^ 

atzz a^zz 

1 



yy '^yz'^zyj 



12 2 I i\ ■ 2 , 
£ — g — e rj \ sin ( 



- e'77'] 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



rjyx — ^xy — (^zz^yx ^yz^zx^ 

f^^zz 

= — i — (—igr/' COS 9b) ■ 
atzz ^ ' 

We write the unit polarization vector as 

^^'^'■^(TT^fT^y^^^^-^'^^-)' 

where iKj = Ex/Ey, iKzj = Ez/Ey, and j is the mode 
index (j = 1 for the X-mode, and 7 = 2 for the 0-mode). 
Eliminating from equation (2.2S), we obtain 



(2.25) 



(2.26) 



1 + 1 



1/2 



(2.27) 



The polarization parameter /3 is given by 

.rjJxx—T]yy 



a = 



e'Vjsin^ 6is + e'V(l-r) 



(2.28) 
2gr]'cos9B 

Using e' — e + a and r;' — i] + a + q [see eq. (2.18)] and 
1 — r — —(m/a) sin^ 9b, we write 



/3 = /3o/3v, 



(2.29) 



where /3o is the polarization parameter in the absence of 
vacuum polarization 



/3o = - 



(e^ — g^ — srj] sin^ 9b 



2grj cos 9b 



1/2 • 2/) / 1 
Mc Sm t/s / 1 + Uc 

I 1 — Iti 



2(1 --uc) cos6»s My 
and /3v is the correction factor due to vacuum polarization 



1 + 



a + q 



1 + 



e{a — q) — rja — erjrn 



(g + m) (1 - Uc) 



1 — lie q+m 



Af-i (1 + % 



For «o 1, equation (2.31) simplifies to 

/3v ~ 



(g + m) (1 - Mo) 



.(2.31) 



(2.32) 



In addition, from equation (2.21), we find 



Kz 



(e — rj — q) sin cos ObKj + g sin 9b 



e sin 9b + (rj + q) cos^ ( 



+ a 



1 — U, 



-—: I sin 6*5 cos 9b Ki 
Ml 



1 - Ui 



— u\^'^Vc sin 9b 



X {(1 - Mc) (1 - u\) (1 + a + gcos^ 9b) 
-Vc [(l-Ui)(l-UcCos^6lfl) -Aftiisin^ 61b] } "\(2.33) 



It is evident from equation (2.33) that the component of 
along k is of order Vc oc p/cu , and thus, at sufficiently 
low densities, Kzj can be neglected so that the modes are 
transverse. 

Finally, the refrac tive index nj of the mode can be ob- 
tained from equation (2.22), which gives 



9V 



e 1 

h -r^ cos 9b 

g K, 



(2.34) 



where tzz ~ e' sin? 9b + rj' cos^ 9b [see eq. (2.20)] and e', rj' , 



g, and a are as given in equation (2. IS 



2.3 Vacuum Resonance 



The polarization parameter /3 [eqs. (2.2J )-(2.32)] directly de- 
termines the characteristics of photon normal modes in the 

medium. For most energies (away from the resonance points, 

1/2 

where /3 — 0; see below), \/3\ cc uj 3> 1, and we have 



r 1 



K2 ^ 2f3 (1/3| » 1). 



(2.35) 



Thus the extraordinary mode [j = 1) is almost linearly po- 
larized with its electric field perpendicular to the k-B plane, 
while the electric field of the ordinary mode (j = 2) is in the 
k-B plane. 

The condition = specifies the resonance points. For 
\f3\ <C 1, equation (2.27) yields 



' sig: 



n(/3) [(-l)^ //^ + |/3|] (1/31 «1). 



(2.36) 



Thus \Kj \ « 1 at /3 = 0, and both modes are circularly po- 
larized. In general, at a given density p, there are three crit- 
ical photon energies at which P = 0. The first is located at 
(jj/iLipe ~ l—iupe/iOBe)'^ /[2{q+m)], i.e., close toui = ujbb- We 
will ignore this critical energy since we are considering only 
the ui <^ LJBe regime in this paper. The second critical en- 
ergy is at u = tOBi', this is simply the ion cyclotron resonance 
and is unrelated to vacuum polarization. The thi rd cr itical 
point ("vacuum resonance") is located at [see eqs. (2.31) and 
(p2 
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Figure 1. The ratio f{B) = Ey{B)/E^\b) as a f unct ion of 
B, where Ey(B) is the vacuum resonance energy [eq. (2.36)] and 



E^\b) = EyiB « Sq) [eq. (|2.3E|)]. 



or at energy 
iiv 



1 a + m 
'^l^'^ 2 



(2.37) 

(2.38) 

m| < 1, 



Note that for 6 < Sn/a p {B < 5 x 10^^ G), |a|, , , 
so that equations ( 2.37) and ( 2.38[) are good approximations 
For B <^ Bq, we use equations ( |2.12[ )-(2.14) to obtain 



Ev{B < Bq) 



E. 



(0) 



rricC 



1.02Fc^/Vi'^Bi~4^ keV. 



(2.39) 



For B ^ Bp, we use the leading terms in equations (2.1C) 
and (2.17) to obtain 

1/2 / R,, X 1/2 



£;v(B > Bq) 



huJr, 



/Stt 

0.69ye^''^Pi/^Bi4'/^ keV. 



(2.40) 

For convenience, we define the dimensionless function f{B) 



via Ev = E^\f{B), so that 



3(5v 



1/2 



(2.41) 



Figure ^ shows that f{B) is a slowly-varying function of b 
lf{B) = 1 for 6 < 1 and f{B) {h/^Y'^ « 0.673 B^f for 
fe > 1; f{B) varies from « 1 at Bm = 1 to 6.7 at Bn = 100]. 

Since By depends on density, a photon with a given 
energy E traveling in an inhomogeneous medium encounters 
the vacuum resonance (/3 = 0) at the density 

pv = 0.964 F-'(Bi4£i)V(B)"'g cm-^ (2.42) 



where E\ — E/{1 keV). In Section 2.4, we discuss how the 
photon polarization state changes as the photon traverses 
the resonant density. 
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Figure 2. The polarization ellipticities [upper pan el: se e 



eq. ( |2.43| )] and refractive indices nj [lower panel; see eq. ( |2.34| )] 
of the photon modes as functions of density near the vacuum 
resonance for B = 5 X G, 0fl = 45°, £ = 5 keV, and Ye = 1. 



2.4 Resonant Conversion of Photon Modes 

In Section 



2.2 



described the photon modes as the 
extraordinary mode (X-mode) and ordinary mode (O- 
mode). This standard way of classifying photon modes (e.g., 
Meszaros 1992) is useful when |/3| ^ 1: the X-mode has 
\Kj\ <^ 1, and its E is mostly perpendicular to the h-B 
plane; the O-mode has \Kj \ ^ 1 and is polarized along the 
h-B plane. The advantage of such a classification is that 
the X-mode and O-mode interact very differently with mat- 
ter: the O-mode opacity is largely unaffected by the mag- 
netic field, while the X-mode opacity is significantly reduced 
(by a factor of order lJ^ from the zero-field value (see 
Section |^. However, the X and O-mode classification be- 
comes ambiguousnear |/3| = i\Kj\ = 1). Indeed, equa- 
tion ( ^.27[ ) or ( ^.36| ) shows that Kj is discontinuous (and 
changes sign) when crossing through /? = at the resonance 
density [eq. (2.42)]. This discontinuity in the mode polar- 



ization vector can be avoided by adopting a diff erent mode 
classification scheme: instead of equation ( 2.27 ), we define 
the plus-mode and minus-mode according to 



1/2 



(2.43) 



Clearly, K± are continuous functions of (3; the plus- mode 
always has K+ > 0, and the minus-mode always has K- < 
0. 

The relationship between the two schemes of mode clas- 
sification (j = 1,2 for the X and O-modes, and J = ± 
for the plus and minus-modes) is illustrated by Figure ^ 
which shows the values of K± and the refractive indices 
n± — ck± /uj for the plus and minus-modes as a func- 
tion of density near the vacuum resonance for E = 5 keV, 
B = 5 X 10" G, and Ob = 45° (see Fig. 1 in Paper II for the 
B = 10^'* G case). The plus-mode (minus-mode) manifests 
as the O-mode (X-mode) at high densities but becomes the 
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X-mode (0-mode) at low densities. If the density variation 
is sufficiently gentle, an O-mode photon created at high den- 
sities will remain on the Jf+'trajectory as it travels outward 
and will adiabatically convert into the X-mode after travers- 
ing the resonance density.^ This is analogous to the MSW 
mechanism for neutrino oscillation (e.g., Bahcall 1989; Hax- 
ton 1995). As shown in Paper II (see also Gnedin et al. 1978; 
Pavlov & Gnedin 1984) ,Q the adiabaticity condition is 



7rcs — 

where 

Sad = 2.52 [/(B) tan 6ls]2/3 1 1 



.|2/3 



1 cm 



(2.44) 

1/3 

keV (2.45) 



and Hp = \dz/d\np\ is the density scale height along the 
ray (evaluated at p = pv). For an ionized hydrogen at- 
mosphere, Hp ~ 2kT /{rripg coaO) = 1.65 r6/((7i4 cos 6) cm, 
where T = lO^Tg K is the temperature, g — 10^''(?i4 cm s~^ 
is the gravitational acceleration, and 9 is the angle between 
the ray and the surface normal. The probability of a non- 
adiabatic "jump" (e.g., from the K+-curve to the K--cmve) 
at the resonance is given approximately by the Landau- 
Zener formula Pjump = exp(— 7r7rGs/2). Thus, for E>2E^d, 
resonant conversion between the X-mode and O-mode is es- 
sentially complete; for E < 0.5i5ad, a, photon will jump across 
the adiabatic curves (Fig. ^), and an X-mode (O-mode) 
photon will remain an X-mode (O-mode) photon in pass- 
ing through the resonance. 

Note that the adiabatic, vacuum polarization-induced, 
resonant mode conversion (from X to O-mode or from O to 
X-mode) is an intrinsically coherent phenomenon: the dis- 
tance over which the conversion takes place is much smaller 
than the photon mean-free path due to absorption or scat- 
tering (Paper II). Such mode conversion is clearly different 
from the incoherent mode-switching due to scattering [see 
eq. (3.3)], which has been included in previous works (e.g., 
Pavlov et al. 1995; Ozel 2001; Paper I). 

As mentioned in Section the formalism of radiative 
transfer in strong magnetic fields developed so far is in- 
adequate for coping with partial mode conversion at the 
vacuum resonance. Thus, in this paper, we shall solve the 
transport equations in the two limiting cases: no mode con- 
version [when the adiabatic condition given by eq. (2.44) 
is never satisfied] and complete mode conversion (when the 
adiabatic condition is always satisfied). In practice, these 
correspond to two different ways of labeling the modes, i.e., 
X, O-modes with j = 1, 2 and plus and minus-modes with 
J = ± (see Sections |^ and ^) 

^ In other words, as a photon propagates across the vacuum res- 
onance under the adiabatic condition, we use the term mode 
conversion to imply that the character of the mode, such as the 
opacity, changes significantly. For example, as a high opacity O- 
mode photon propagates from high to low density, it becomes a 
low opacity X-mode photon. Alternatively, this adiabatic regime 
could be called mode conservation since the photon remains on 
the same plus-mode (or minus-mode) branch across the reso- 
nance. We use the former terminology. 

^ In the context of radio wave propagation in inhomogeneous 
plasmas, the analogous problem of linear transformation of modes 
has been extensively studied (e.g., Budden 1961; Ginzbur g 197 0V 
^ Bezchastnov et al. (1996) apparently used equation (2.43) as 
the basis of the mode; therefore they automatically included 



2.5 Mode Collapse and Breakdovirn of Faraday 
Depolarization 

Here we comment on the phenomenon of mode collapse 
(see also Gnedin & Pavlov 1974; Pavlov, Shibanov, & 
Yakovlev 1980; Soffel et al. 1983; Bulik & Pavlov 1996), 
which appears under more res trict ive conditions than the 
reso nance s discussed in Section 2.3 [see discussion following 
eq. ( 2.36 )]. Mode collapse occurs when the two polariza- 
tion modes beco me id entical, i.e., Ki = K2 or K+ = K- 
[eqs. (2.27) and (2.43)], which requires P — iiy'r ~ ±j. 
Obviously this is only possible when the dissipative terms, 
i.e., those involving As = 1 -f ius/uj, in the dielectric t enso r 
[eq. (ETI)] are retained. Also, we see from equation (2.34) 



that the two modes have the same indices of refraction when 
mode collapse occurs. 

To de termine the mode collapse point, we use equa- 
tion (^.28|), i.e.. 
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- — e' rj (1 + m/a) sin^ 



(2.46) 



2gr\' cos 9b 

Note that e' ,g, r( , and fi are considered complex in this sub- 
section. The condition Re(/3) = yields the three critical 
energies E'^ ~ E\, Ebe, Ebz, and these are the resonant en- 



ergies discussed in Section 2.3. Mode collapse a lso re quires 
Im(/3) = ±1, so that (3 = ±i. From equation (2.46), it is 
easily shown that 



Im(/3) = /37e 



cos 9f 



(2.47) 



where /3 depends on the photon energy E — hui, magnetic 
field B, and density p and 7e = Ve/ui = 7c, -I- 7r is the 
dimensionless damping rate of the electron. For a hydrogen 
plasma, the collisional and radiative damping rates are 



^ 9.24 X W-''piT~^^^E^^ (1 



e-^/''^) (2.48) 



7r 



9.52 X lO'^-Ei, 



(2.49) 



respectively, where we have neglected the Gaunt factor. The 
damping rate of the ion is 71 « 7cme/mp. Clearly, at a given 
E (« Ey , Ebb, Ebi), the condition Im(/3) = ±1 selects out 
a critical angle 6b = ^coii at which mode collapse occurs. 
For example, for p = 1 g cm"^ and B = 10^" G, \i3\ ^ 
350,1.3 X 10*^,0.21 aX E EY,EBe,EBi, respectively, and 
for B = 10^^ G, ]/3] « 5.6,130,0.63.0 Since |/3|7o < 1 at 
E ~ Ev and i?si for these parameters, we see that the 
mode collapse associated with the vacuum resonance and 
the ion cyclotron resonance occur at 9b = 6'coii close to 90° . 

More generally, the modal description of radiative trans- 
fer is valid only in the limit of Faraday depolarization, i.e., 
when the condition |Re(n+ — n_)| ^ |Im(n-|- + n_)| is sat- 
isfied (Gnedin & Pavlov 1974). One can show that away 
from the mode collapse points, this is indeed the case. The 
breakdown of Faraday depolarization at _E = EBi appears 
to be concentrated near 9b = 90° , but the situation for the 
vacuum resonance is more complicated. For Mo 3> 1 (and 
neglecting the ion effect), we find 



(complete) mode conversion. Ozel (2001) used equation (2.27) 



as the basis and therefore did not include the mode conversion 
effect. 



At E Ey, one can show that (3 ~ — u 
in the limit Ue ^ 1 and Ey ^ EBi- 



1/2 



/2 = -EBji2Ey) 
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Re(n+ 

Im(n+ 
where 



-1/2 



COS ^sRe 



1/2 .2, 

Uc sm ( 



2 cos9p 



9s + ito ( 1 + cos 



0] 



- «7cj 



(2.50) 
(2.51) 

(2.52) 



and /3v ~ 1 - (g + = 1 - {E / EyY [see eq. (|2.3l| )] 



Thus, for l/9v|^7c (i.e., for i5 slightly away from -Ev), 
Faraday depolarization is satisfied for all Ob- On the 
other hand, for |/?v| 7c (i-e., for E extremely close 
to -Ev), the condition |Re(n+ — n_)|>jlm(n+ + n_)| re- 
quires (7ciiy^/2)| sin^ 6b/ cos6b\ ~ llm(/3)| < 1; this condi- 
tion translates to 9b ^ ^coii or 6s > vr — Scoii- In other words, 
the breakdown of Faraday depolarization at the vacuum res- 
onance (|/3v| ^ 7e) is restricted to an angular region of 
Ob around 7r/2. To satisfy Faraday depolarization for most 
Ob requires 7eUe'''^/2^''^ ^ 1, where 7e and Ue are evalu- 
ated at the vacuum resonance. For a given magnetic field 
B, this requirement implies E ^ i5Far(B): photons with 
energy E much greater than E-pry{B) encounter the vac- 
uum resonance [at py(E)\ at which Faraday depolarization 
is satisfied for most 6b\ conversely, for E < i5Far(-B), Fara- 
day depolarization breaks down at the vacuum resonance 
for a significant range 6b around n/2. Figure ^ shows 
EFar(B) as a function of B. Note that when the ion ef fect is 
i nclud ed, the simple expressions given by equations (2.5C)- 
(2.52) are modified, and we evaluate the complex n+ and 



using equation (2.34) and obtain E-pz.^{B) numerically 



[for i3i4>30, the function Ep^r{B) is multi- valued because 
Re(n^ — n?.) depends non-monotonically on E]. For exam- 
ple, at B = 5 X 10" G, photons with £;<£Far ^ 2 keV 
encounter the vacuum resonance between the X-mode and 
O-mode decoupling layers [po < Pv{E) < p^^^ see Sec- 
tion ^ and Fig. ^ , and Faraday depolarization breaks down 
at this resonance for most 6b- 

Currently it is not clear how significant an effect the 
breakdown of Faraday depolarization will have on the ra- 
diative transfer. Given that this breakdown of the modal 
description occurs for such a narrow range of E [i.e., /3v ^ 
7o ^ 1; note that this range is even narrower than the width 
of the opacity feature around Ey, see Section |^ and eq. (3.1t ) 
in particular], the effect should be small; but one cannot be 
sure. The above analysis reinforces our statement in Sec- 
tions 1^ and 2.4 that a rigorous treatment of radiative trans- 
fer near the vacuum resonance, in general, requires solving 
the transport equations for the four Stokes parameters. We 
plan to study this and other related issues in a future pa- 
per. We shall ignore the mode collapse phenomenon and the 
breakdown of Faraday depolarization in our calculations in 
the rest of the paper. 
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> 



T = 5xl0'5 K 
0„ = 45° 



with ions^ 



0.1 



10" 



lOifs 
3 (G) 



lO'" 



Figure 3. Critical energy _EFar(S) (at T = 5 X 10^ K and 
6b = 45°) above which the condition for Faraday depolariza- 
tion is satisfied at the vacuum resonance p = o\r( E). The solid 



lines indicate BFar calculated from equation (2.34) with ions in- 
cluded, while the dashed lin e indi cates a pprox imate values of -Epar 
calculated from equations (2.5C) and ( |2.5l| ) which neglect ions. 
Note that at large magnetic fields (_B > 3 X 10^^ G), the function 
SFar(B) (with ion effects included) is multi- valued. The unshaded 
region is where Faraday depolarization is valid, while the shaded 
regi ons are where Faraday depolarization breaks down (see Sec- 
tion for more details). 



[1 ± {Kj cos 6b + Kz,j sin 6b) 



2{l + K]+Kl^) 



sm Ob 



i cos Ub ) 



1 + 



(3.1) 



(3.2) 



wher e Kj is given by equation (2.27) for j = 1,2 or equa- 
tion (2.43) for j — ±, and Kz,j is given by equation (2.33). 

The electron scattering opacity from mode j (= 1,2 for 
the X and O-modes or ± for the plus and minus modes) 
into mode i is given by (Ventura 1979; Ventura, Nagel, & 
Meszaros 1979) 



1 

p ^ (UJ + aUJBey + Vl 
a — — 1 



(3.3) 



where ric is the electron number density, ctt is the Thomson 
cross-section, and yl^ is the angle integral given by A^^ — 
(3/87r) J dfe |e^j^. The electron scattering opacity from mode 
j is 



3 BEHAVIOR OF OPACITIES AND VACUUM 
POLARIZATION RESONANCE 

The radiative opacities depend on the normal mode polar- 
ization vector through its projection on the rotating frame 
with the z-axis along B. The cyclic components of e-* are 



V2 



p ^-^ {io -)- aiOBe 

a — — 1 



(3.4) 



where A^ = X/i-i ^q- transverse-mode approxima- 

tion {Kz,j = 0), the polarization vector e? satisfies the 
completeness relation X/^-i l^±l^ ~ (l -I- cos^ 6b) 12 and 
X^j^i jfizl'^ ~ sin'^Ss, and thus Aa = 1. For bj ^ tJse, a 
suppression factor {ijj/ujb&)'^ in the opacity results from the 
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strong confinement of electrons perpendicular to the mag- 
netic field. Similar features appear in the electron free-free 
absorption opacity (e.g., Virtamo & Jauho 1975; Pavlov & 
Panov 1976; Nagel & Ventura 1983) 



ff,o _(^sr^ ^ I J ,2 ff 

J' ^ n ^ '■• • ■- ^2 ' •■2 1 "I 



(3.5) 



where 



Z CL^ — 2- -t:^ — TV^^ ) 



ff3\/3 1 



. -ff ' (3-6) 

rii is the ion number density, and af axe the free-free 

absorption coefficient and velocity-averaged free-free Gaunt 
factor, respectively, in the absence of a magnetic field. In 
equation (p^sj), g±i = g^ and gg = are the velocity- 
averaged free-free Gaunt factors in a magnetic field, which 
we evaluate using the expressions given in Meszaros (1992) 
(see also Nagel 1980). 

The ion contribution to the scattering and absorption 
opacities are 



'7-2 \ ^ 
Z rric \ nifjT 

Arrir, 



(3.7) 



Z rric \ niCTT 



E; 



p (t^ — O-LOBiY + 

a— — 1 



^CK i (3.( 



1 / Z^rue 



2 1 
QO \ ^ 



(c- 



+ vt 



I J |2-ff 



(3.9) 



when crossing the resonance), the opacities are determined 
with j = ± (right panels) .n Note that since iJv depends 
on p, the left and right panels of Fig. ^ represent genuinely 
different opacities in the radiative transfer (see below). 

The behavior of the opacities near the ion cyclotron 
resonance Esi = 0.63 B\a keV (for Z — A = 1) was studied 
in Paper I. The vacuum resonance feature at Ev [eq. ( 2.40 )] 
arises from the interplay between the p lasm a and vacuum 
polarization effects discussed in Section 2.3 (e.g., Meszaros 
& Ventura 1979; Pavlov & Shibanov 1979; Ventura, Nagel, 
& Meszaros 1979) . To understand the behavior of the X and 
O-mode {j = 1,2) opacities near the vacuum resonance, we 
write the scattering and absorption opacities as 



■n-oO-T 



with 



(3.10) 



+ 



(1-^)^ + 7? 



J l2 



(1 + v^r 



,(3.11) 



where M' = {Amp/Zmc)Z~'^^^ for scattering and M' = 
{Amp/Zma)Z^^^ for absorption, 71 = Vi/to, and we have set 
the free-free Gaunt factors to unity in Kf°^ and Aa = 1 in 
Kfj^ for simplicity. Away from the ion cyclotron resonance, 
the ion contribution to the opacities is negligible, and ~ 
(|e^+P + |ei|^)/-uc + Hl^- In the limit of > 1 (and using 
the transverse a pprox i mati on Kz. i ~ or Ve <^ 1), we find 



from equations (2.35), (3.1), and (3.2) that 



1 |2 



1 =F - cos 6b 



(2/3)- 



■sin Ob (3.12) 



Note that the ion cyclotron resonance occurs for a = +1, 
i.e., when the electric field of the mode rotates in the same 
direction as the ion gyration. The total scattering and ab- 
sorption opacities in a fully ionized medium is then the sum 
of the electron and ion components, namely 
and nf" = ^.f" + 



In equations (3.3)-(?.9), we have included damping 
through fc — fr,o + fc.e aud = u^.i + i/c,i, where ;/r,c ~ 



(2e /3meC )a; and Ur. 



(Z ma/Amp)iyr,c are radiative 



damping rates and i^cc = (aosf /'T'cO"T)fr,c and Uc,i — 
(rric/ Amp)Vc,e are collisional damping rates (see Pavlov et 
al. 1995, and references therein). For the photon frequencies 
of interest, u) 3> Vc^i, and therefore damping is negligible 
except near resonance. 

Figure ^ shows examples of the absorption opacity nf'^ 
as a function of photon energy for mag netic field B = 10^'' G 
and Ob = 45° at densities and a temperature characteristic 
of NS atmospheres. For angles 6b not too close to 0° or 180° 
(e.g., 10° <6'b < 170°), the opacity nf^OB = 45°) is indica- 
tive of the behavior of K|'"^(fe), while K,f^^{k) exhibits strong 
angle dependence for 6b near 0° or 180°. Because of r nod e 



conversion due to the vacuum resonance (see Section 2.4), 
two sets of photon modes are depicted in Fig. ^ when mode 
conversion is neglected (e.g., an X-mode photon remains in 
the X-mode when traversing the vacuum resonance), we use 
j = 1,2 for the X and O-modes (left panels); on the other 
hand, if we assume that mode conversion is complete for all 
energies (e.g., a plus-mode photon remains in the plus-mode 



I 2 i2 1 2 a 
\C'±\ = -^COS dB 



1 ± 



/3 cos I 



■ 2 a 

sm 6 b ■ 



(3.13) 



Thus, for the X and O-modes,^ 

~ — + 7^9 sin'^^s, ^-sin'^s (l/3| > 1). (3.14) 

On the o ther h and, in the limit of ^ 1, we find from 
equation (2.36) that 



J |2 



J\2 



i{[l±(-l)^sign(/3), 



(-1)^1/91 + TT 



sin dB 



sin^ 6b- 



(3.15) 
(3.16) 



^ In plot ting the ri ght pa nels of Fig. ^, we have used /3 = |/3o|/9v 
[see eqs. (2.29) and (2.30)] in evaluating the polarization param- 



eter f3 so that the plus or minus-mode is resonant at the ion 
cyclotron energy. If we use (3 = /3o/3vi then the -|- and — curves 
will switch across EBi since /3 changes sign. However, since Ebi 
is independent of p, it is easy to show that the radiative transfer 
is the same for both cases (/3 = \Po\Pv versus [3 = /3o/3v)- 
^ Note that in this section, we are considering 9b away from 
and it; when 6b is close to or tt, the dominant terms in are 
different. 
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.1 1 10 0.1 1 10 

E (keV) E (keV) 

Figure 4. Absorption opacity k^^^ as a function of energy for B = 10" G, Bb = 45°, and T = 7 X 10^ K at p = 370g cm ^ (upper 
panels) and p = 0.1 g cm~^ (lower panels). The left panels show the opacities for the X-mode (j = 1, solid lines) and the O-mode {j = 2, 
dashed lines), while the right panels are for the plus- mode {j = +, solid lines) and the minus-mode (i = — , dashed lines). In the left 
panels, the dotted lines show the the X-mode opacity when vacuum polarization is neglected. Several features are marked in the figure: 
is the ion cyclotron energy, Ey is the vacuum res onan ce energy, and Epe is the electron plasma energy. The variations of the O-mode 



opacity at i? ^ Epe are due to K^j [oc ujpf,; see eq. ( 2.33)] 



Thus we have 

singes [i + {-iyw\] 



1/31 « 1). 



(3.17) 



Clearly the X-mode {j — 1) opacity is significantly enhanced 
near the vacuum resonance (/3 ~ 0), while the O-mode opac- 
ity is only reduced by a factor of two from its usual value. 
At the resonance, the two modes become circularly polarized 
and possess the same opacities (see left panels of Fig. ^) . 
Equations (3.14) and (3.17) determine the "line shape" 



of the vacu um r es onan ce feature. For x = (E — E\)/ Ev <^ 1, 
equations (E^2m)-( 2.32) imply |/9v| ~ 2la;| and 



with 



l/2|i 
Ur 1- 



, sin^ 6b 
cos 9 b 



cos Ob 



u' 1- 



us \ sm 



2 , 



(3.18) 



(3.19) 



Thus 



i sin^ Qb (1 — la::l/a::v) 
u^^ + (a;v/2x)^ sin^ Qi 



\x\ < Xy 

Xy <C l^l <C 1 



.(3.20) 



This should be compared to the value when there is no vac- 
uum polarization (nv), 



^(nv) 



-1 , 2 . 2 

-I- xv sm Ob 



(3.21) 



The strength of the vacuum resonance feature can be mea- 
© 0000 RAS, MNRAS 000, 1-|^ 



sured by the dimensionless "equivalent width" 
1 



Tv 



n/j^ 1-jy J res 



KxdE 



SX •^■•cs J -1/2 



1/2 



ixdx, 



(3.22) 



where the integral is over the region in which Cx/d""' > 1, 
i.e., Ey/2 <E< 3Ey/2 or -1/2 < a; < 1/2 (see Fig. ^. Using 
equation (3.20) and for sin^ 6b ~ cos6b ~ 1, we find 



,1/2 



1 - Ui 



(3.23) 



It is important to note that the region \x\ < ( a few xy) <^ 1 
contributes most to the integral in equation ( |j.22| ). This is 
illustrated by the shaded region in Figure |^, which shows a 
schematic pictu re of the vacuum resonance feature; note that 
from equation ( 3.20 ), the continuum is uZ^ lower than the 
peak. This region must be resolved in numerical calculations 
of atmosphere spectra (see Section^. Also note that at low 
densities where Ey lies below the ion cyclotron energy so 
that Uj > 1, the vacuum resonant feature is narrower [see 
eq. ( |3.1£| ); also compare the two left panels in Fig. |^. 

When mode conversion is assumed to be complete, the 
plus and minus-mode opacities k± do not exhibit a line fea- 
ture but rather show a step function-like change at the vac- 
uum resonance energy Ey (see the right panels in Fig. M). 
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x = (E-Ey)/E^, 



0.5 



Figure 5. A schematic picture of the X-modo opacity near 
the vacuum resonance energy Ey. The dashed Une represents 
the opacity when vacuum polarizatio n is neglected. The shaded 
region, —3xy ^ x ^ Sxy [see eq. (3.19); note that xy <^ 1], 
provides the dominant contributio n to the e quiva lent width of 
the vacuum resonance feature [eqs. (3.22) and ( |3.23| )], although a 
much broader region, —0.5 ^ x ^ 0.5, is affected by vacuum po- 
larization. The opacity at the resonance {E = Ey) is ^ Uc times 
larger than the continuum. 



4 QUALITATIVE DISCUSSION OF THE 
VACUUM POLARIZATION EFFECT ON 
ATMOSPHERE EMISSION 

In order to understand qualitatively the effects of vacuum 
polarization on the radiation spectra from magnetar atmo- 
spheres, we estimated in Paper II the location of the decou- 
pling layer at which the optical depth is of order unity for 
photons of different energies and polarization modes. Here 
we conduct a somewhat more accurate calculation by in- 
cluding the energy dependence of the Gaunt factor and show 
results for diflFerent magnetic field strengths. For simplicity, 
we consider fully ionized hydrogen atmospheres, adopt the 
transverse approximation {K^j — 0), and present results 
for only 9b ~ 45°. We also treat the temperature as a con- 
stant, which is a good approximation when estimating the 
decoupling density since T varies at most by a factor of a 
few while p varies by many orders of magnitude above the 
decoupling layer. For g = {GM/R^){1 - 2GM/Rc^)-^^^ = 
2.4 X 10" cm s"^ (corresponding to a M = 1.4 M©, R = 
10 km NS), hydrostatic equilibrium and the ideal gas equa- 
tion of state, P — 2pkT/mp, give the density p at column 
density y (in g cm~^) as p = 1.45 Tg"^?/ g cm~'^, where 
Ts = r/(10^ K). The optical depth (= / K^dy) at p is 
then given by Tv = 0.69 Te Kv dp' . Note that we use here 

the effective opacity defined by = \Kf"'{i<^'' + ^i'J')]^'^ 
(see Rybicki & Lightman 1979). The decoupling density p 
can then solved for by setting the optical depth = 1. 
Numerical results of the decoupling density are shown in 
Figure | for B = lO" and 5 x lO" G. 

We can obtain analytic expressions to understand the 




E (keV) 

Figure 6. The density at which photons of different polarization 
modes decouple from the matter as a function of photon energy 
for T = 5 X 10" K, Bg = 45°, and B = IQi* G (lower panel) and 
B = 5 X lO^** G (upper panel). The solid lines show the X-mode 
decoupling density p^"^' when vacuum polarization is included 
but mode conversion is neglected, the short-dashed lines are for 
p^^' with no vacuum polarization, the long-dashed lines are for 
the O-mode po (note that po is unaffected by the magnetic field 
and vacuum polarization effects), and the d ot-dashed lines show 
the vacuum resonance density py [eq. (2.45)]. The vertical dotted 



lines deno te the critical energy E^^ for adiabatic mode conversion 
[eq. ( ^.45| )]. In general, vacuum polarization (including the mode 
conversion effect) reduces the decoupling density of X-mode pho- 
tons with energies between min(i?ci , i?ad) an d E r9.. where E^i 
and Ec2 are given approximately by equations ( [4.4| ) and ( [4.3| ), re- 
spectively. The ion cyclotron resonance Egi = 0.63 Bn is clearly 
evident, though it can be suppressed when vacuum polarization 
is taken into account, as in the case of iJ = 5 X 10^* G. 



various features in Fig. ^ by setting the Gaunt factor to 
unity and by ignoring scattering, which is only important 
for > 10 keV, even though the results shown in Fig. ^ do 
account for the actual variation of the Gaunt factor and 
use the effective opacity which includes scattering. First, 
let us ignore vacuum polarization and proton effects. The 
absorption opacity can then be written as Kf°^ « kqCj > where 



-1/2 



E'-^S cm^g-^ [with 5 = 1 - e 



-E/kT 



Pi = 



Ko ~ 9.3piTg 
g cm" 

and Co ~ 1 and ~ u^^ [see eqs. (3.10) and (3.14) 



p/(l g cm ■^), and Ei = E/{1 keV)] is the zero-field opacity 



The decoupling densities of the two modes [with no vacuum 
polarization (nv)] are then approximately 



(nv 

Po 

(nv 

Px 



0.56 Tfi 



-1/4^3/2 



-1/2 -3 

' g cm 



1/2 (nv) 
uj p'q ' . 



(4.1) 

(4.2) 
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Thus the X-mode photons emerge from deeper in the atmo- 
sphere and are the main carriers of the X-ray flux. 

Next, we include vacuum polarization and protons but 
neglect resonant mode conversion. For the O-mode, the 
opacity is largely unaffected by the magnetic field and vac- 
uum polarization effect. Thus the decoupling density po is 



still /5q . For the X-mode, when p^^' > pv or 



77 t: 



1/6^-2/3 



/(B)^/^ keV 



(4.3) 



the photons created at p'^^^^ will encounter the vacuum res- 
onance, near which the X-mode opacity is greatly enhanced; 
thus the decoupling density [with no mode conversion (nc)] 



P^x'^^ will be smaller than Px^'- For E>Ec2, p''x^' is close 
to p^x^\ with the small difference due to the opacity be- 
ing slightly modifled by the vacuum effect even away from 
the resonance (see the left panels of Fig. To evaluate 
the X-mode optical depth At across the resonance region 



(nc) 



1/2 K 

(e.g., 0.9 < p/pv < 1.1) we note that /3 ~ u^Jx [eq. ( 3.18 )1 
for |a::| = |Ap/pv|^0.1. Using equations (3.14) and (3.17), 

we find ~ l/{uax'^) for Ue^^'^ ^ |a;|<0.1 and ^jc ~ 1 

x\ <C Uc^^"^ [compare with eq. ( 3.20 )]. We then have 



for 



Ar ~ 0.7 fZoTapvu-^'^ = QATI'^pIE- 



'1/2 



S, which is 



1/2 

a factor ~ larger than the optical depth of the non- 

resonant region (/9<0.9/9v)- Thus at energies below iJc2, 
p^x'^^ closely follows pY until E drops below another critical 
energy Ec\ , which is set by Ar ~ 1 or 

£ci ~ 10Tg-'/*B-/''V(B)' keV. (4.4) 

Below Eel, Px'^^ returns to approximately p^^^ (see Fig. |^. 

Now consider the effect of mode conve rsion at the vac- 
uum resonance. For E > 1.3 -Ead [see eq. ( 2.45 )], adiabatic 



mode conversion is nearly complete (Pjump < 3%). The O- 
mode photons traveling from high densities through py are 
converted to X-mode photons, which then freely stream out 
of the atmosphere. Thus, for 1.3i5ad < Ec2, the effective 
decouphng density for the X-mode is px = pv- For E > Ec2, 
the vacuum resonance occurs inside both the X and O-mode 
decoupling layers, and therefore px ~ Px^"*- For E <^ -Ead, 
mode conversion is ineffective; therefore px = Px'^'' ■ Around 
-Ead, the X-mode photons are emitted from both p^'^' (with 
probability Pjump) and pv [with probability (1 — Pjump)]. 
Also note that Px'^^ ^ pv for E > Ed . 

Figure ^ and the above analysis show that vacuum po- 
larization reduces the decoupling density for photons with 
energies between min(Eci , -Ead) and -Bc2, i.e., these photons 
decouple in shallower, lower temperature regions of the at- 
mosphere. This gives rise to a rather broad depression fea- 
ture in the emergent radiation between min(-Eci, -Ead) and 
-Ec2 (see also Paper II). In other words, even though the vac- 
uum resonance opacity feature is very narrow, it produces a 
broad depression because of its density dependence and the 
large density gradient present in NS atmospheres. Clearly, 
to quantify the depth of the depression, one must solve for 
the temperature profile of the atmosphere self-consistently 
(see Sections |^ and ^ . 

Finally, the sharp absorption feature in Fig. ^ is due to 
the ion cyclotron resonance at Esi ~ 0.63 B14 keV. Previ- 
ous atmosphere models which neglect vacuum polarization 
yield relatively large equivalent widths for the proton cy- 



clotron line (Paper I; Zane et al. 2001). However, we see 
from Fig. M that when Esi lies within the depression trough, 
i.e., Ed <EB^ < £c2 or 3 X lO" G < BfiBy^^T^^^'' < 2 x 
10^^ G, the continuum fiux around Ebi is greatly reduced,^ 
and the equivalent width of the ion cyclotron line becomes 
much narrower. This suppression of the ion cyclotron line by 
vacuum pol ariza tion makes the line more difficult to detect 
(see Section 6.2). 

In Appendix ^ we present a toy atmosphere model 
that mimics the vacuum resonance effect discussed in this 
section. 



5 NUMERICAL METHOD FOR TREATING 
VACUUM POLARIZATION 

We solve the full, angle-dependent radiative transfer equa- 
tions for the two coupled photon modes in order to con- 
struct self-consistent NS atmosphere models. Because the 
method for rigorously treating partial mode conversion has 
not been developed, we consider two limiting cases: no mode 
conversion and complete mode conversion at all energies 
(see Section 2.4). In our models, the temperature corrections 
AT(r) at each Thomson depth r are applied iteratively un- 
til AT(r)/r(T) < 1%, deviations from radiative equilibrium 
are < 1%, and deviations from constant ffux are < 2% (see 
Paper I for details of our numerical method). Here we dis- 
cuss the numerical difficulties involved when attempting to 
include the effects of vacuum polarization in the atmosphere 
models. 



5.1 No Mode Conversion: Grid Resolution 

As shown in Section ^, vacuum polarization induces a nar- 
row, density-dependent resonance feature in the X-mode 
opacities (see Figs. ^ and|^. The "equivalent width" of the 
vacuum resonance feature at -Ev(p) is dominated by the nar- 
row energy range AE/Ey — \E — Ey\/Ey <?> xy ^ uZ^^'^ ^ 
10-^Ei/Pi4, where Ei = E/{1 keV). In a real atmosphere, 
density and energy gradients are smooth so that such a nar- 
row and density-dependent feature is resolved and accounted 
for completely. Numerical atmosphere models, however, nec- 
essarily require finite grids of discrete depth (density), en- 
ergy, and angle. The difficulty with a finite grid resolution 
is illustrated in Figure ^. Consider a photon with energy E 
that happens to lie very close [|-E — -Ev(pi)l/E^v(pi) ^v] 
to the vacuum resonance energy Ev{pi) of one of the den- 
sity grid points pi; such a photon will encounter a region of 
strongly enhanced opacity and will therefore decouple from 
the matter at a shallow depth. On the other hand, a pho- 
ton with energy E', which lies outside the resonance regions 
of all density grid points will not encounter any enhanced 
opacity regions and will decouple deeper in the atmosphere. 



' If we use E^^<Egi, then the first inequality becomes 4 X 
IQl-* G <B/(B)-2/3(tan6»fl)-2/3. Note that the range of mag- 
netic fields must be considered very approximate since Ed, i?ad, 
and Ec2 depend on the direction of photon propagation, and re- 
distribution of photon spectral fiux occurs in real atmospheres. 
For example, even for B = 10^'* G, the width of the ion cyclotron 
line is reduced by vacuum polarization (see Section ^.2| ). 
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Figure 7. A schematic picture of the radiative opacity as a 
function of energy at three neighboring density grid points Pi_i, 
pi, Pi -1-1. Density and depth increase upwards, and the vacuum 
reson ance feat ure o ccurs at a higher energy in deeper levels [see 
eqs. ( 2.39 ) or (2.40)]. The opacity at a given depth has been offset 
for clarity. In the "equal grid" method, the energy grid points 
En-i, En, -En+i are chosen to be equal to the vacuum resonance 
energy at one of the density grid points, i.e.. En = Evipi)- The 
long-dashed line for pi indicates the apparent opacity obtained 
by inter pola ting energy grid points E„—i, En, and E„+i (see 
Section p.l\ ). The energy E lies close to one of the resonance 
energies, while E' lies outside any resonance region. 



Therefore, if we choose a general energy grid that is not tied 
to the density grid, the computed radiation spectrum will 
exhibit narrow absorption lines at energies which happen to 
lie close to Ev at one of the density grid points. Such nar- 
row lines are numerical artifacts of the finite grid resolution 
since, as discussed in Section ^, vacuum polarization is ex- 
pected to only produce a broad depression in the emission 
spectrum. 

Clearly, to correctly account for the density- 
dependent vacuum resonance, one must have a depth 
grid that is sufficiently dense so that the opac- 
ity features at neighboring grid points overlap ap- 



preciably, i.e., [Ev{pt+i) - Ev{pi)]/Ey{pi)<xv or 

{pi+i — pi)/pi<Uc^^'^ ~ 10~'^EiB^_^. Such a dense depth 
grid is beyond the capabilities of the current generation of 
atmosphere models. 

One method to avoid artificial "absorption lines" in the 
computed spectrum is to adopt an "equal grid" method: 
we use an equal number of depth (density) and energy 
grid points with every energy grid point being placed at 
En = Ev{pi) (see Fig. ^). The advantage of this method is 
that it guarantees a smooth spectrum. However, this method 
has problems as well. In order to determine the atmosphere 
structure through temperature corrections (see Paper I, and 
references therein), it is necessary to compute mean opaci- 
ties, which involve integrations over energy, at a given depth; 
using the "equal grid" method leads to overestimates of the 



strength of the opacity feature (see Fig. ^ the area under 
the long-dashed line overestimates the area under the solid 
line, which represents the true opacity). Nevertheless, by in- 
creasing the number of depth grid points, the "equal grid" 
method will yield atmosphere models increasingly close to 
reality (see Section ^) . 

An alternative method to produce smooth spectra in- 
volves "saturating" the opacities near the vacuum resonance 
(Oz el 2 001). We shall comment on this method in Sec- 



tion 6.3 



5.2 Complete Mode Conversion 

If we ass ume the adiabatic condition given by equa- 
tion ( |2.4^ ) is satisfied for all energies, then a plus-mode 
(minus-mode) photon will remain in the plus-mode (minus- 
mode) as it traverses the vacuum resonance (see Section 2.4). 
In this "complete mode conversion" limit, we simply solve 
the coupled radiative transfer equations for the plus and 
minus-mode photons rather than the X and 0-mode pho- 
tons. The numerical procedure is analogous to the one out- 
lined inPaper I. The modes are calculated using equa- 
tion ( |2.43[) for the polarization parameter rather than equa- 
tion ( p.27| ). In this case, the vacuum resonance manifests 
as a step function-like feature in the opacities, as discussed 
in Section ^ and shown in the right panels of Fig. ^, and 
there is no numerical difficulty in handling such an opacity 
feature. 



6 NUMERICAL RESULTS 
6.1 Atmosphere Structure 

The temperature profiles for fully ionized hydrogen atmo- 
spheres with Toff = 5 X 10^ K and magnetic field B — 
5 X 10^'' G oriented perpendicular to the surface (i.e., the 
angle between the field and the surface normal is Ob = 0) 
are plotted in Figure ^. The curves marked "nc#" are mod- 
els which include vacuum polarization but no mode con- 
version and using the equal grid method described in Sec- 
tion 5.1 , with the 



"#" indicating the number of grid points 
per decade in Thomson depth r. When vacuum polarization 
is neglected (nv), the temperature profiles show a plateau at 
r ~ 1 — 100 (such a feature was already noted in Paper I^ 
and in Ozel 2001^ ). This arises because the X-mode pho- 
tons decouple at r ~ a few x 10^ and thus have very little 
energy exchange with the matter, while for t>1, the O- 
mode photons are ineffective at regulating the temperature 
because of their small ffux. We see from Fig. ^ that vacuum 
polarization tends to diminish this plateau in the temper- 
ature profile since vacuum polarization causes the X and 
O-mode opacities to be similar to each other (see Fig. ^ 
and the decoupling densities to be closer (see Fig. ^). We 



° The "no vacuum polarization" models presented here differ 
slightly from those in Paper I because of a numerical error in our 
calculation of the Gaunt factors in Paper I. This error has been 
corrected in the present paper. 

® Ozel's models include vacuum polarization; see Section 6.3. 
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Figure 8. Temperature T as a function of Thomson depth r for 
fully ionized hydrogen atmospheres with i? = 5 X 10^"* G, ©s = 0, 
and T(,ff = 5 X 10® K. The curves marked "nc#" are models which 
include vacuum polarization but no mo de conversion using the 



"equal grid" method described in Section 5.1, with "#" indicating 
the number of grid points per decade in t. The long-dashed line 
is for t he sa me model as "nc20" but restricting /3v > 10~^ (see 
Section |6.3| ), while the dot-dashed line is for the model with no 
vacuum polarization (nv). 



also see that vacuum polarization tends to increase the tem- 
perature at most depths (and the general temperature gra- 
dient) relative to the case without vacuum polarization. In 
the deepest layers where both modes are diffusing, a rise in 
the X-mode opacity due to the vacuum resonance feature 
reduces the bandwidth over which flux can be transported, 
and a larger temperature gradient is required to maintain a 
constant total flux; these layers are therefore backwarmed 
to higher temperatures (see Mihalas 1978 for a general dis- 
cussion of backwarming) . On the other hand, in the plateau 
r ~ 1 — 100, the free streaming X-mode photons, generated 
from deeper, hotter layers, encounter a sharp rise in opacity 
and are absorbed, which causes the temperature to increase. 

Figure ^ compares the temperature versus density pro- 
flies for various magnetic fields and vacuum polarization ef- 
fects. For the models which include vacuum polarization but 
neglect mode conversion (nc) , we use twenty grid points per 
decade in r. Compared to models without vacuum polariza- 
tion (nv), the models with complete mode conversion (mc) 
show higher temperatures in the deep layers: although the 
O-mode photons at a particular depth are converting into 
X-mode photons and carrying away heat, X-mode photons 
from deeper, hotter layers are converting into O-mode pho- 
tons and depositing more energy. We also note that larger 
magnetic fields cause the vacuum resonance to occur at a 
lower energy for a given dens it y or a t a hi gher density for 
a given energy [see eqs. (2.39), ( |2.40| ), and (2.42)] and thus 
shift the region that is heated by vacuum polarization effects 
to deeper layers. 

In Appendix we construct toy atmosphere models 
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Figure 9. Temperature T as a function of density p for fully 
ionized hydrogen atmospheres with T^g = 5 X 10® K, 0^ = 0, and 
B = 10" G (lower panel), B = 5x 10^* G (middle panel), B = 
10^® G (upper panel). The solid lines are for models with vacuum 
polarization but no mode conversion (nc), the dashed lines are 
for models which assume complete mode conversion (mc), and 
the dot-dashed lines are for models with no vacuum polarization 
(nv). 



with opacities which mimic the vacuum polarization effect. 
The temperature profiles of the toy models are qualitatively 
similar to the results discussed in this section. 



6.2 Spectra 

Figure |l^ shows the spectra of fully ionized hydrogen at- 
mospheres with B — 5 X 10^* G, 0_g = (magnetic field 
perpendicular to stellar surface) , and Tcs = 5 x 10® K. Plot- 
ted alongside is the blackbody spectrum at T = 5 x 10® K. 
This figure illustrates the effect of grid resolution for models 
which include vacuum polarization but no mode conversion. 
As discussed in Section 



5.1 



a lower grid resolution tends to 
overestimate the strength of the vacuum resonance feature, 
which leads to the stronger depression at ~ 0.2 — 2 keV. In 
the following, we shall use the "nc20" model when we show 
results which include vacuum polarization but no mode con- 
version and adopt the abbreviation "nc". 

Figures h_]J-ll3 depict the spectra for fully ionized hy- 
drogen atmosphere models with B = 10^"* G, 5 x 10^* G, 
and 10^^ G, respectively, Ob = 0, and Teff = 5 x 10® K. In 
each figure, the spectra from three models are depicted to- 
gether with a blackbody spectrum at T = 5 x 10® K: a model 
which neglects vacuum polarization (nv), a model which in- 
cludes vacuum polarization but no mode conversion (nc), 
and a model which assumes complete vacuum polarization- 
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Figure 10. Spectra of fully ionized hydrogen atmospheres with 
B = 5 X IQi-* G, 013 = 0, and Tefi = 5 X lO'^ K. The curves 
marked "nc#" are models which include vacuum polarization but 
no mod e co nversion using the "equal grid" method described in 
Section 



5.1 



vith indicating the number of grid points per 

decade in r. The corresponding temperature profiles are shown 
in Fig. |. The dotted line is for a blackbody with T = 5 X 10'^ K. 
Note that the absence of the ion cyclotron feature at Ebi for 
model nc6 is an artifact of the low grid resolution. 



induced mode conversion (mc). Recall that the nc model 
and the mc model represent two limiting cases, and the true 
results with partial mode conversion (see Section 2.4) are 
expected to lie between the nc and mc curves. We see from 
Figs. pT| - [l^ that the differences in the spectra between the 
two limiting cases are not significant. 

Figures |l^-|l3| show that, when vacuum polarization is 
neglected (nv), the spectra exhibit significantly harder high 
energy tails and a depletion of low energy photons relative to 
the blackbody. The hard tails were already noted in previous 



studies of NS atmospheres with _B < 10 G (e.g., Shibanov 
et al. 1992; Pavlov et al. 1995; see also references cited in 
Paper I), and they arise because high energy photons have 
lower opacities and thus decouple from deeper, hotter lay- 
ers. We see from Figs. pl[]l3| that vacuum polarization can 
significantly reduce these high energy tails (£ > a few keV) 
and causes the spectra to be softer compared to the mod- 
els without vacuum polarization. This reduction of the high 
energy tail is due to the broad depression caused by the 
density-dependent vacuum resonance feature, as discussed 
in Section^ (see Fig. 

Figures h_]J-ll3 also reveal that vacuum polarization can 
significantly suppress the ion cyclotron line feature in the 
spectra. When vacuum polarization is neglected, the ion cy- 
clotron line is broad (see Paper I for more discussion of the 
feature; see also Zane et al. 2001), and the equivalent widths 
are 0.14 keV {Esi = 0.63 keV), 2 keV (Ssi = 3.15 keV), and 
4 keV [Ebx ~ 6.3 keV) for B14 = 1, 5, and 10, respectively. 
When vacuum polarization is included but mode conversion 
is neglected (model nc), the equivalent widths become 20 eV 
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Figure 11. Spectra of fully ionized hydrogen atmospheres with 
B = IQl-* G, 0B = 0, and Teff = 5 X 10^ K. The solid line is for 
an atmosphere with vacuum polarization but no mode conversion 
(nc), the dashed line is for an atmosphere with complete mode 
conversion (mc), the dot-dashed line is for an atmosphere with no 
vacuum polarization (nv), and the dotted line is for a blackbody 
with T = 5 X 10'^ K. The Esi = 0.63 keV ion cyclotron feature 
(from ^ 0.4 to ~ 1 keV) has an equivalent width ~ 140 eV 
when vacuum polarization is not included (model nv) but is only 
20 eV in the no mode conversion case (model nc) and ~ 6 eV 
in the complete mode conversion case (model mc). 



{Esi = 0.63 keV), 0.1 keV {Eei = 3.15 keV), and 0.2 keV 
{Ebz = 6.3 keV), i.e., the equivalent widths of the ion cy- 
clotron line have been reduced by a factor of ~ 10; when 
mode conversion is assumed to be complete (model mc), the 
equivalent widths are even smaller: 6 eV {Ebz = 0.63 keV), 

0. 09 keV {Esi = 3.15 keV), and 0.03 keV (Ssi = 6.3 keV), 

1. e., a reduction by a factor > 10. This reduction of the 
ion cyclotron line strength by vacuum polarization was ex- 
pected from the analysis in Section ^ (see Fig. and it 
occurs when the ion cyclotron energy Ebz overlaps with the 
broad depression caused by the vacuum resonance. The re- 
duced width of the ion cyclotron line makes the line difficult 
to observe with current X-ray detectors (see Section]^. 

A comparison of the spectra for the models with differ- 
ent magnetic fields is shown in Figure ^ for the case which 
includes vacuum polarization but no mode conversion (nc) 
and in Figure ha for the case which includes vacuum po- 
larization with complete mode conversion (mc). Note that 
despite the suppression of the high energy tails by the vac- 
uum polarization effect, all spectra are still harder than the 
blackbody spectrum. 



6.3 Comparison with Previous Work 

As mentioned in Section |l|, there have been few previous 
works on magnetar atmosphere models that include the ef- 
fect of vacuum polarization. The most recent and relevant 
one is that by Ozel (2001), who constructed self-consistent 
atmosphere models with parameters similar to those de- 
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Figure 12. Spectra of fully ionized hydrogen atmospheres with 
S = 5 X IQi" G, 0ij = 0, and T^g = 5x 10*^ K. The solid line is for 
an atmosphere with vacuum polarization but no mode conversion 
(nc), the dashed line is for an atmosphere with complete mode 
conversion (mc) , the dot-dashed line is for an atmosphere with no 
vacuum polarization (nv), and the dotted line is for a blackbody 
with T = 5 X lO^ K. The Esi = 3.15 keV ion cyclotron feature 
(from ~ 2 to ~ 5 keV) has an equivalent width ~ 2 keV when 
vacuum polarization is not included (model nv) but is only ~ 
0.1 keV in the no mode conversion case (model nc) and ^ 0.09 keV 
in the complete mode conversion case (model mc). 
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Figure 13. Spectra of fully ionized hydrogen atmospheres with 
B = 1015 G, 0B = 0, and Teff = 5 X 10^ K. The solid line is for 
an atmosphere with vacuum polarization but no mode conversion 
(nc), the dashed line is for an atmosphere with complete mode 
conversion (mc), the dot-dashed line is for an atmosphere with no 
vacuum polarization (nv), and the dotted line is for a blackbody 
with T = 5 X 10^ K. The Esi = 6.3 keV ion cyclotron feature 
(from 5 to ~ 10 keV) has an equivalent width ^ 4 keV when 
vacuum polarization is not included (model nv) but is only ^ 
0.2 keV in the no mode conversion case (model nc) and 0.03 keV 
in the complete mode conversion case (model mc). 



scribed in our paper, i.e., T^a ~ 5x 10*^ K, B ~ lO^^-lO^^ G, 
Ob = 0, and the atmosphere consists of fully ionized hydro- 
gen. While there are qualitative resemblances between some 
aspects of our results and Ozel's, e.g., the temperature pro- 
files show plateau features which are weakened by vacuum 
polarization (compare our Fig. |^ to her Fig. 5), there are 
also major differences in the physical ingredients and meth- 
ods used in our models and those of Ozel's. (1) Ozel ap- 
plied the vacuum polarization formulae which are valid only 
for B<Bq = 4.4 X 10" G. At B > 10^^ G, these incor- 
rect formulae underestimate the vacuum resonance energy 
by a factor of a few (see Fig. 0). (2) Ozel neglected ions 
in the plasma response and the ion cyclotron resonance in 
the opacities. As discussed in Paper I, the ion effects can- 
not be neglected since they influence the spectrum in the 
same energy range as vacuum polarization effects. Figure 
compares the spectra of atmosphere models which include 
ions and neglects ions. Despite the suppression of the ion 
cycl otro n feature by the vacuum polarization effect (Sec- 
tion |6.2[ ) , the "with ion" and "no ion" cases have apprecia- 
bly different spectra (e.g., at _B = 10 keV, the flux of model 
"nc" differs from model "nc/no ions" by a factor of about 
ten). (3) Ozel's work neglects the mode conversion effect 
due to vacuum polarization (Paper II; Section 2.4). (4) Even 
when mode conversion is neglected, Ozel's method of treat- 
ing the va cuum resonance differs from ours. As discussed 
in Section 




5.1 



to fully account for the density-dependent 
vacuum resonance feature in the opacity (and thus produc- 
ing smooth spectra) , one would need to have a prohibitively 



E (keV) 

Figure 14. Spectra of fully ionized hydrogen atmospheres with 
0B = E'lid TeS = 5 X 10^ K when vacuum polarization is 
included but mode conversion is neglected (nc) . The long-dashed 
line is for the B = 10^'^ G atmosphere, the solid line is for the 
B = 5 X 10^* G atmosphere, the short-dashed line is for the 
B = 10^'' G atmosphere, and the dotted line is for a blackbody 
with T = 5 X 10'^ K. 
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Figure 15. Spectra of fully ionized hydrogen atmospheres with 
©B = and Teff = 5 X 10^ K when mode conversion is complete 
(mc). The long-dashed line is for the B = 10^^ G atmosphere, the 
solid line is for the i? = 5 X 10^* G atmosphere, the short-dashed 
line is for the B = 10^^ G atmosphere, and the dotted line is for 
a blackbody with T = 5 X 10*' K. 



high density /depth grid resolution. Ozel produces smooth 
spectra by adopting a "saturation" method, whereby the 
peak X-r node opacit y is restricted by requiring j/3vj > 0.01 
[see eqs. (2.32) and (3.18); recall that the vacuum resonance 
occurs at /3v = 0]. However, as discussed in Section ^, the 



II < 



-1/2 



0.001 at 



region with |/3vj ~ \E/Ew 
1 keV for B = lO" G; see eq. ( ^^ ) and Fig. || contributes 
most to the "equivalent width" of the vacuum resonance 
feature. Therefore the saturation method does not capture 
the main contribution of the vacuum resonance feature in 
the opacities and significantly underestimates the effect of 
vacuum polarization. Figure pi shows the temperature pro- 
file of the atmosphere model using our "equal grid" method 
together with the restriction j/3v| > 0.01. Clearly, since the 
saturation scheme does not capture the entire effect of the 
opacity change, the heating of the atmosphere (compared to 
the case in which vacuum polarization is neglected) due to 
the vacuum resonance feature is less than the model with 
no saturation restriction. Figure ^ compares the spectra 
with and without saturation. Two effects are responsible for 
the qualitative differences. First, because saturation under- 
estimates the "equivalent width" of the vacuum resonance 
feature, it does not yield the true decrease in flux (see Sec- 
tion ^) so that the the broad depression does not reach as 
low of energies. Second, the softer spectra of the satura- 
tion model at energies E^l keV is the result of the lower 
temperature at the photon decoupling layers in this model 
(see Fig. ^). We also note again that even though vacuum 
polarization produces harder spectra at high energies than 
the blackbody spectrum at the same effective temperature, 
the spectra are still softer than atmosphere models which 
neglect vacuum polarization (see Figs. |ll|-[l3[). 

It is appropriate to reiterate that our own numerical 



Figure 16. Spectra of fully ionized hydrogen atmospheres with 
_B = 5 X IQl-* G, 013 = 0, and T^ff = 5 X 10** K. The solid line is for 
an atmosphere with vacuum polarization but no mode conversion 
(model nc, which includes ions), the dot-dashed line is for an 
atmosphere with vacuum polarization but no mode conversion 
and neglecting ion effects (nc/no ions), the long-dashed line is 
for an atmosphere with no vacuum polarization (model nv, which 
includes ions), the short-dashed line is for an atmosphere with no 
vacuum polarization and neglecting ion effects (nv/no ions), and 
the dotted line is for a blackbody with T = 5 X 10'^ K. 
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Figure 17. Spectra of fully ionized hydrogen atmospheres with 
B = 5 X IQl* G, Gs = 0, and Teg = 5 x lO'' K. The soUd line 
(no saturate) is for an atmosphere with vacuum polarization but 
no mod e co nversion using the "equal grid" method described in 
Section 5.1, the dashed line (saturate) is for the same atmosphere 
with vacuum polarization but no mode con version using the "sat- 
uration" method described in Section |3.3[ and the dotted line is 
for a blackbody with T = 5 x 10^ K. 
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treatment of the vacuum polarization effect also has limi- 
tations, and these are discussed in Section ^ We also note 
that Ozel's calculations include the full angular dependence 
of the scattering opacity, whereas the results presented in 
this paper are based upon an approximate calculation of 
the scattering source function in order to speed up the com- 
putation (see Section 2.2 of Paper I). Our numerical tests 
show that this approximation does not affect the temper- 
ature profile of the atmosphere and only produces a small 
reduction (by less than a factor of two) in the ffux at very 
high energies {E^ 10 keV). Finally, we have not examined 
the angular dependence and beaming pattern of the emission 
as done in Ozel (2001) (see also Ozel, Psaltis, & Kaspi 2001 
and Ozel 2002 for analyses of the resulting pulse profiles and 
pulse fractions). 



7 DISCUSSION AND CONCLUSION 

We have presented a detailed study on the atmospheres 
and thermal radiation of isolated neutron stars with super- 
strong magnetic fields B > 10^* G and effective temper- 
atures Toff ~ a few X 10^ K. Such a study is needed since 
surface emission has already been detected from a number of 
magnetar candidates (AXPs and SGRs), and current/future 
observations, when combined with theoretical modeling, can 
potentially provide important constraints on the properties 
of these enigmatic objects (see Section ^ . 

Following up our previous work on magnetized neutron 
star atmospheres (Paper I), we focus on the effect of vac- 
uum polarization in this paper. It was already known from 
earlier theoretical studies (see Section |l| and Meszaros 1992 
for references), which we have generalized to the B > 10^** G 
regime and included the effect of ions, that vacuum polariza- 
tion changes the dielectric property of the plasma and gives 
rise to a resonance feature in the opacity at photon energy 
Ev « 1.02Y^^^^pY^B-lf{B) keV, where = Z/A is the 
electron fraction, pi = p/(l g cm"'^), B14 = 5/(10^* G), 
and f{B) is a slowly-varying function of B of order unity 
(see Fig. |l|). Furthermore, it was shown in Paper II (see 
also Gnedin et al. 1978; Pavlov & Gnedin 1984) that pho- 
tons with energies E> a few keV propagating in the at- 
mospheric plasma can adiabatically convert from one polar- 
ization mode into another at the vacuum resonance den- 
sity pv X 0.96Y~^EfBf4f{B)-^ g cm-^ where Ei = 
E/{1 keV); this resonant mode conversion greatly influences 
the radiative transfer because the two modes have vastly dif- 
ferent opacities in strong magnetic fields. In this paper, we 
have attempted to incorporate these vacuum polarization ef- 
fects into self-consistent atmosphere models. Through both 
analytic considerations (Section ^; see also Paper II) and nu- 
merical calculations (Section P), we have shown that vacuum 
polarization leads to a broad depression in the high-energy 
(i?> a few keV) radiation fiux from the atmospheres as 
compared to models without vacuum polarization. Despite 
the rather sharp vacuum resonance feature in the opacity 
at a given density, the depression in the spectrum is broad 
(from a few keV to tens of keV, depending on the field 
strength) because of the large density gradient in the at- 
mosphere. As a result, the atmosphere emissions possess 



softer high energy tails than models without vacuum po- 
larizationrj. 

Another important effect of vacuum polarization on the 
atmospheric spectra is that the strength of the ion cyclotron 
line is greatly suppressed when vacuum polarization is in- 
cluded in the atmosphere models (see Figs. |ll[ - |l^ ). This 
is a direct consequence of the aforementioned flux depres- 
sion caused by vacuum polarization (see Fig. when the 
depression trough overlaps with the ion cyclotron line at 
EBi = 0.63ye-Bi4 keV, i.e., for 10^* G<B<5 x 10^^ G, the 
continuum flux around Esi is greatly reduced, and the ion 
cyclotron line appears narrower (see Section W. For exam- 
ple, we flnd that the equivalent width of the 3.15 keV proton 
cyclotron line for B = 5 x lO^** G is reduced from about 
2 keV (no vacuum polarization) to about 0.1 keV (with vac- 
uum polarization). Obviously the reduction of the ion cy- 
clotron line width makes the line more difficult to detect. 
Indeed, recent observations of several AXPs with Chandra 
and XMM-Newton X-ray telescopes failed to resolve any sig- 
nificant line features in the spectra (e^., Patel et al. 2001; 
Juett et al. 2002; Tiengo et al. 2002)[J this could possibly 
be indicating the presence of vacuum polarization effects. 

We note that the radiative transfer formalism adopted 
in this paper (and in previous atmosphere models by other 
researchers) relies on the transfer of two photon polariza- 
tion modes. This is inadequate for treating the vacuum 
polarization-induced resonant mode conversion effect, espe- 
cially because the effectiveness of m ode conversion depends 
on photon energy (see Section 2.4 and Paper II). The at- 



mosphere models studied in this paper represent two limit- 
ing solutions to the transfer problem (no mode conversion 
and complete mode conversion), and they are expected to 
bracket the true solution. Nevertheless, to properly account 
for the mode conversion effect associated with vacuum polar- 
ization, as well as mode collaps e an d the breakdown of Fara- 
day depolarization (see Section 2.5), one must go beyond the 
modal description of the radiation field by formulating and 
solving the transfer equation in terms of the photon density 
matrix (or Stokes parameters) and including the effect of a 
nontrivial refractive index. We plan to study this problem 
in the future (see also Lai & Ho 2002b). 

We caution that our models, as well as previous models 
of magnetar atmospheres (see Section ^), assume that the 
atmospheres are completely ionized. For the magnetic field 

For all the models studied here, the spectra are still harder 
than the blackbody spectrum (with the same effective tempera- 
ture) because of the non-grey opacities. Also note that in this pa- 
per we are only concerned with thermal emission from the neutron 
star surface. Nonthermal emission or the reprocessing of the ther- 
mal emission by the magnetospheric plasma may give rise to high 
energy, power-law tails in the magnetar spectra (e.g., Thompson, 
Lyutikov, & Kulkarni 2002). 

" For AXP 4U0142-I-61, Juett et al. (2002) give an upper limit of 
about 10-50 eV for the equivalent width of any line feature in the 
spectrum; this is comparable to our predicted value of ~ 6 — 20 eV 
at B = IQi* G and ~ 100 eV at B = 5 x 10^** G. In addition, 
it should be noted that the spectra presented in our paper corre- 
spond to emission from a local patch of the neutron star surface. 
When the spectra are integrated over the entire observable sur- 
face, which necessarily involves (very uncertain) variations in the 
magnetic field strength and direction and effective temperature, 
the equivalent width of the cyclotron line will be further reduced. 
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strengths and surface temperatures characteristic of AXPs 
and SGRs, it is not clear that bound atoms and molecules 
have sufficiently small abundances to contribute negligibly 
to the opacity (see the discussion section in Paper I), given 
their greatly enhanced binding energies in strong magnetic 
fields (see Lai 2001 for a review). However, it may be the 
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case thiit, like the ion cyclotron hue, absorption Imes due to 
bound apecies are Huppreayed If they fall within the broad 
flux depression caused by vacuum polarization. The absence 
of any resolved features in recent observations of several 
AXPs (Patel et al. 2001; Juett et al. 2002; Tiengo et al. 2002) 
may therefore be indicative of the vacuum polarization effect 
in superstrong magnetic fields. Also of concern is the dense 
plasma effect on the radiative transfer (see Paper I) . Clearly, 
much work remains to be done before we can have complete 
confidence in our quantitative understanding of the radia- 
tion from the surfaces of strongly magnetized neutron stars. 
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APPENDIX A: ATMOSPHERE INCLUDING 
VACUUM POLARIZATION RESONANCE 
EFFECTS: TOY MODELS 

Al Opacities 

In Section ^ we examined the effect of the vacuum polar- 
ization resonance feature on the depth at which the photons 
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Model (a) 



Mod 



j-2 



j = l 



(b) 



Figure Al. Toy model opacities used in Appendix ^ with 
j = 1,2 indicating the X and O-mode, respectively, and J = it in- 
dicating the plus and minus-modes, respectively. Model (a) mim - 
ics the case where mode conversion is neglected [see eq. (^^]; 
Model ( h) m imics the case where mode conversion is complete 
[see eq. ( |A3| )]. The opacities ar c p arameterized by the location of 
the resonance feature i?v [eg. ( |A2| )] , the width parameter A, and 
kh and kl- 



that comprise the surface spectrum are emitted. To deter- 
mine the emission spectrum, the temperature profile must 
be determined self-consistently. In this section, we consider 
several toy atmospher e mo dels based on the simplified opac- 
ities shown in Figure 



Al 



These models serve to illustrate 
the key physical effects of the vacuum resonance feature and 
the numerical subtleties when constructing real atmosphere 
models (Sections ^ and 

As shown in Section |^, in a real atmosphere, the O- 
mode {j = 2) opacity is largely unaffected by the vacuum 
resonance while the X-mode {j — 1) opacity exhibits a sig- 
nificant peak around E\r. Therefore, for Model (a), we as- 
sume the O-mode opacity is grey, k£ — ku, and the X-mode 
opacity has a simple square function form: 



X 
Ke = 



where 

Ev = 



for E < Ev or E > \Ev 
for Ey < E < XEy 



CP 



1/2 



[Model (a)], (Al) 



(A2) 



p is i n g cm ^, E is in keV, ( ^ B^l f{B) [see eqs. (2.39)- 
(2.41)], and the width of the line feature is (A — 1)-Ev- Equa- 



tion (Al) thus resembles the X-mode opacity shown in the 
left panels of Fig. ^ near the vacuum polarization resonance. 
From Section ^ we know that /kl ~ Uc at the resonance 
peak, and the width is (A — 1) ~ {u\^^\l — ?ii|)~^; we will 
see that only the "equivalent width" (kh/'«l)(1 — A~^) ~ 
uj'^/\l - Mil (~10^ at 1 keV for B = lO" G) matters for 
the radiation spectrum. Model (a) mimics the real atmo- 
sphere models where vacuum-induced mode conversion is 
neglected. 




E (keV) 

Figure A2. The upper panel shows the optical depth (in 
terms of the opacity kl) at which photons of different modes 
decouple from the matter as a function of energy for T = 5 X 
10^ K, C = 0.3, A = 1.3, kh = Kg=, and kl = 10-*Kg=. The lower 
panel shows the corresponding density at these depths. The 
solid lines are for the X-mode with an opacity given by Model (a) 
(the O-mode decoupling depth = 1 is not shown), an d the 
short and long-dashed lines are for Model (b) (see Fig. The 
vertical d otte d lines indi cate the critical energies i?ci [eq. ( |A5| )], 
E,2 [eq. (0)], E+ [eq. (O], and E- [eq. (113)]. 



In Model (b), we assume the plus and minus-mode opac- 
ities are simple step functions: 



kh for E < Ey 

kl for E > Ev 

kl for E < Ev 

kh for E > Ev 



[Model (b)] 
[Model (b)]. 



(A3) 



These opacities resemble the behavior of the plus and minus- 
modes in the right panels of Fig. ^. Thus Model (b) mimics 
the atmospheres where mode conversion is assumed to be 
complete. 



A2 Photon-Matter Decoupling Region 

As in Section ^, we can calculate the density pd where pho- 
tons of different modes decouple from the matter. The cor- 
responding decoupling depth Td(-B) is determined from 1 = 
f^'^ dTL[K,-E.{TL) / t^i], where cLtl = —pni^dz. Similar to Sec- 
tion ^ we assume a constant temperature profile when cal- 
culating Td and pd (see Section ^ for justification). Thus, for 
a fully ionized hydrogen atmosphere, p{tl) = 1.45 j/Tg"^ = 
PotlTq^ g cm~^, where po ~ 1.45/kl (here kl is in cm^g~^ 
and po is in g cm~^). Combining the density profile with 
equation (Al) gives the X-mode opacity as a function of tl, 
i.e., KE(rL) = KH for {Te/po){E/XCf < tl < {Te/po)iE/Cf 
and ke(tl) — kl otherwise. We then obtain the decoupling 
depth Td of the X-mode as a function of energy, 
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ME) 



1 - (fCH/ACL - 1) (l - p-iTs {E/Cf 

for E < E^i 
«:l/«:h + (1 - kl/kh) Pq 'Tb [S/ (CA)]^ 

for Eel < E < Ec2 
1 for _B > £c2 



.(A4) 



The two critical energies, -Eci and Ec2, corresponding to 
when the decouphng depth rd — {Ta / pi)){E / C,Y and rd = 
{n/po){E/\Cf are 



1 



-1/2 



i?c2 = CPo''r6"''"A. 

The decoupling density, pd 

poTg-i - (kh/«:l - 1) (l - A- 

for S < E^i 
poTg'^AiL/^H + (1 - kl/kh) [E/ {CX)f 



p(rd), is given by 

■ {E/0^ 



P^iE) = 



(A5) 
(A6) 



(A7) 



for Ed < E < Ec2 
PoTg-i for E > E,2 



Note that when kl/kh ^ 1 and A ~ 1, the decoupling 
density pd(-E) ~ {E/Cf = pv{E) for E^i < E < Ec2, which 
is in agreement with the result of Section ^ (see Fig. |^. 
Equations (A4) and (A7) are plotted in Figure ^2 with ^ = 
0.3, A = 1.3, Kii = k"^" = 0.4 cm^ g-\ and kl = IQ-^ 



the 



choice of these parameter values is discussed in Section A3 
The minimum decoupling depth and density occur at Ed 

1 



Pd{Ec 



A-2+kh/kl(1-A-2) 



PoTg 



(A8) 



(A9) 



A-2+kh/kl (l-A-2)' 
Clearly the minimum decoupling depth depends critically on 
the equivalent width of the opacity feature, i.e ., (kh /kl)(1 — 
A~^), which corresponds to Fv in equation (3.22). In fact, 
in the limit of kl/'^h ^ 1 and A « 1, the parameters kl, 
kh, and A enter equations (A4)-(A£) only through the com- 
bination (kh/kl)(1 — A~^). The larger the equivalent width, 
the shallower the decoupling depth and hence the cooler the 
region where the observable photons are produced and the 
deeper the depression t hat is produced in the emission spec- 
trum. Also, equations (|A5|) and (A6) show that, by increas- 



ing the magnetic field (or decreasing Q, the broad depression 
is shifted to lower energies. 

Now consider Model (b) with opacities given by equa- 
tion (|A3|). Using the same procedure as above, we obtain 
the decoupling densities for the plus and minus-mode: 



pZ(s) = 

where 

Et = 
E- = 



Po 



T^^f^L/fiH + (1 - kh/kl) (E/Cf 



for E < EZ 
for E > E+ 



PoTs' - {Ku/tiL - 1) (E/(f 

for E < E^ 
PoT^'^i^l/kh for E > E~ 



^ l/2rr~l/2/ I n1/2 

CPo (fi:L/fi:H) . 



(AlO) 



(All) 



(A12) 
(A13) 



Equations (AlO) and (All), along with the corresponding 
decoupling de pths for the plus and minus-modes, are also 
plotted in Fig. A2 . 



A3 Numerical Method 

Using the opacities of Models (a) and (b) in our magnetic 
diffusion atmosphere code (for simplicity and computation 
speed) described in Paper I, we obtain self-consistent atmo- 
sphere models for a given choice of A, kh, and kl. The 
results presented in this section are somewhat modified for a 
full, angle-dependent solution of the radiative transfer equa- 
tion, but since here we are more interested in the qualitative 
effects of vacuum polarization, the diffusion models are ade- 
quate for illustrative purposes. Temperature corrections and 
deviations from radiative equilibrium and constant fiux are 
<1% (see Paper I). 

For illustrative purposes, we choose the parameters 
A, «H, and kl for the toy models that roughly correspond 
to the case of an atmosphere with magnet ic fie l d B = 5 x 



10 G. For this magnetic field, equations (|2.39|)-(2.41) give 



C ~ Biif(B) « 0.3. We choose a fairly broad width (A - 
1) = 0.3 for the feature so that the required grid resolution is 
not too computationally demanding. We also choose values 
for the parameters kh ~ t^o and kl = 10~*Ko'' so that the 
equivalent width of the feature {k.h/k,-l){1 — X~^) — 4000, 
which corresponds roughly to the equivalent width of th e 
vacuum resonance feature for B = 5x 10^* G [see eq. (3.23)]. 

As discussed in Section |^, to correctly account for the 
density-dependent opacity feature, such as in Model (a), in 
numerical calculations, it is important that the density or 
depth grid is sufficiently dense so that the opacity features 
at neighboring density grid points overlap. For Model (a) 
with A = 1.3, the minimum required grid spacing is | Ap|/p < 
(A- 1)2 = 0.09. In our calculations, we use 50 grid points per 
decade in r (the total number of depth grid points is then 
D ~ 350 covering seven decades in r), which corresponds 
to (pi — pi^\)/pi ~ 0.04 — 0.05. In addition, to obtain an 
accurate representation of the opacity feature in the energy 
domain, we choose an energy grid that is determined from 
the depth grid: for each depth grid point, we place an energy 
grid point on either side of the opacity edge, so that the 
total number of energy grid points A*' = 4_D for Model (a) 
and TV = 2D for Model (b). In this way, we completely 
characterize the shape and strength of the opacity features 
in Models (a) and (b). 



A4 Numerical Results 

Figure plots the temperature profiles for the atmosphere 
models with = 5 x 10^ K using Model (a) (mimicking 
the case of no mode conversion) and Model (b) (mimicking 
the case of complete mode conversion) and C, — 0.3, A = 1.3, 
Kh = kq'*, and kl ~ 10~*kq''. Also shown is the tempera- 
ture profile from a double-grey atmosphere, i.e., two photon 
modes each with constant opacities Kg = kl for the X-mode 
and k£ = Kh for the 0-mode. We see that both Model (a) 
and Model (b) have higher temperatures (at all depths) than 
the double-grey case. This is due to the heating effect asso- 
ciated with the increase in opacity (see Section g and Figs. ^ 
and^. Also note the temperature inversions that occur near 
T ~ 600; this is due to the absorption of high-energy X-ray 
photons as they encounter the resonance feature. Figure A4 
shows the corresponding spectra of the atmosphere models. 
Note that all of these spectra are harder than a blackbody 
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Figure A3. Temperature T as a function of Thomson deptii t for 
atmosph ere models using the toy opacities given in Appendix ^ 
(see Fig. 0) and T^ff = 5 X 10^ K, C = 0.3, A = 1.3, kh = 
and kl = IO^^Kq"^. The solid line is for the temperature profile of 
an atmosphere with Model (a), the long-dashed line is the profile 
for Model (b), and the short-dashed line is for a double grey 
atmosphere. 

at T = 5 X 10*^ K and are similar to the behavior seen in 
Fig. FL3, except without the ion cyclotron line. 
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Figure A4. Spectra of atmospheres with T^ff = 5 X 10® K, 
C = 0.3, A = 1.3, Kii = Kq', and kl = IO^^Kq*^. The solid line is 
for an atmosphere using the Model (a) opacities, the long-dashed 
line for an atmosphere using the Model (b) opacities, the short- 
dashed line is for a double grey atmosphere, and the dotted line 
is for a blackbody with T = 5 X 10** K. 
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